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ABSTRACT 


Incoherent  averaging  of  radar  returns  measured  at  different  frequen¬ 
cies  has  been  shown  to  be  an  effective  method  of  improving  fixed  target 
detection  performance.  This  improvement  results  from  the  fact  that  frequen¬ 
cy-averaging  narrows  the  amplitude  probability  distributions  of  both  the 
target  and  clutter  returns  (after  averaging),  thereby  allowing  returns  from 
^  targets  to  the  more  easily  distinguished  from  those  of  clutter. 

This  report  examines  one  aspect  of  this  improvement  -  the  effect  of 
i  frequency  averaging  on  decreasing  the  sample  variance  (i.e.,  narrowing  the 

distribution)  of  the  measured  clutter  returns.  The  attainable  reduction  in 
clutter  sample  variance  that  results  from  frequency-averaging  is  shown  to  be 
strongly  dependent  upon  the  spatial  "non-unformity”  of  the  clutter,  as 
characterized  by  the  variance  of  its  intrinsic  spatial  amplitude  distribu¬ 
tion.  The  more  non-uniform  the  clutter,  the  less  effective  frequency-aver¬ 
aging  will  be  in  reducing  the  clutter  sample  variance  and  corresponding  CFAR 
threshold  setting.  An  analytic  expression  is  derived  for  the  amount  of 
sample  variance  reduction  provided  by  frequency  averaging,  as  a  function  of 
clutter  spatial  standard  deviation  and  the  number  of  frequencies  averaged. 
The  amount  of  variance  reduction  predicted  by  this  expression  is  comfirmed 
*  by  experimental  measurements  (at  Ku  band)  of  clutter  returns,  at  16  and  32 

different  frequencies,  for  a  variety  of  different  clutter  types.  The  sepa¬ 
ration  between  the  frequencies  at  which  these  measurements  were  made  was  on 
the  order  of  the  radar  bandwidth,  thereby  ensuring  that  the  returns  measured 
at  those  frequencies  were  statistically  independent  of  one  another. 


iii 


For  nose  commonly  encountered  types  of  ground  clutter,  it  is  found 
that  averaging  clutter  returns  over  four  different  frequencies  reduces  the 
clutter  sample  standard  deviation,  and  the  corresponding  CFAR  threshold  off** 
set,  by  35%  to  50%  relative  to  that  required  for  data  measured  at  a  single 
frequency  (no  averaging).  However,  little  additional  reduction  in  threshold 
level  appears  to  be  obtained  when  the  returns  are  averaged  over  more  than  8 
frequencies. 


iv 


CONTENTS 


Abstract 

List  of  Illustrations 
List  of  Tables 
1 .0  INTRODUCTION 

1.1  Computation  of  CFAR  Detection  Threshold 

1.2  Effect  of  Frequency-Averaging 

1.3  Organization  of  Report 

2.0  THE  TWO-PARAMETER  CFAR  ALGORITHM 
3.0  EVALUATION  OF  FALSE  ALARM  PROBABILITY 

4.0  EFFECT  OF  FREQUENCY-AVERAGING  ON  THE  CFAR  THRESHOLD  SETTING 

4.1  Effect  of  Logarithmic  Detection 

4.2  Clutter  Sample  Variance  as  a  Function  of  Frequencies 
Averaged  and  Intrinsic  Clutter  Variance 

4.3  Effect  of  Frequency-Averaging  on  Threshold  Setting 
5.0  Experimental  Verification 

5.1  The  Variance  Reduction  Factor 

5.2  Description  of  Experimental  Clutter  Measurement 
Data 

5.3  Estimation  of  Intrinsic  Spatial  Clutter  Standard 
Deviation  From  Measured  Clutter  Sample  Statistics 

5.4  Comparison  of  Experimental  Data  with  Theoretical 
Predictions 

ACKNOWLEDGMENTS 


REFERENCES 


Appendix  A  Derivation  of  CFAR  Stencil  Variance  as  a  Function 
of  Stencil  Size,  Number  of  Frequencies  Averaged, 
and  Clutter  Statistics 

Appendix  B  Normality  of  the  Probability  Distribution  of  the 
Frequency  Averages  of  Logarithmically  Detected 
Clutter  Returns 

Appendix  C  Variance  of  the  Output  of  a  Logarithmic  Detector 
in  Response  to  a  Rayleigh  Distributed  Noise  Input 


LIST  OF  ILLUSTRATIONS 


Diagram  of  a  typical  2-dimensional  CFAR  stencil, 
illustrating  the  threshold  detection  process 

Probability  distributions  of  target  and  clutter  returns 
for  a  given  value  of  target-to-clutter  ratio,  indicating 
effect  of  detection  threshold  setting  on  target  detection 
and  clutter  false  alarm  probabilities.  Effect  of 
frequency-averaging  on  narrowing  of  the  distributions 
is  also  illustrated. 

Probability  distribution  of  log-detected  clutter  amplitudes 
with  and  without  frequency-averaging  illustrating  the 
effect  of  frequency-averaging  on  narrowing  the  distribution 
and  reducing  the  threshold  level  required  to  maintain  a 
specified  Pfa. 

Effect  of  frequency-averaging  in  reducing  clutter  sample 
standard  deviation  (square-root  of  sample  variance  of 
frequency-averages  over  space)  for  different  clutter 


LIST  OF  TABLES 


5.1 

Summary  of  ClutCer  Parameters*,  Stockbridge,  NH  - 

28 

Late  Summer,  1977 

5.2 

Summary  of  Clutter  Parameters*,  Ft.  Devens,  MA  - 

29 

Summer,  1978 

5.3 

37 

H 


« 


* 


viii 


1.0 


INTRODUCTION 


The  detection  of  fixed  ground  targets  in  clutter  is  based  upon  the 
amplitude  of  the  return  from  a  target  exceeding  that  from  the  surrounding 
ground  clutter.  The  detection  process  used  is  that  of  the  CFAR  (Constant 
False  Alarm  Rate)  detector,  in  which  the  amplitude  of  each  resolution  cell 
is  compared  against  a  threshold,  This  threshold  is  set  at  a  specified  level 
above  the  average  amplitude  of  the  clutter  in  the  vicinity  of  the  resolu¬ 
tion  cell  being  tested  for  the  presence  of  a  target  return.  If  the  ampli¬ 
tude  of  the  return  in  the  resolution  cell  being  tested  exceeds  this  thresh¬ 
old,  a  target  is  declared  to  exist  in  that  cell. 

However  the  amplitude  of  the  clutter  returns  is  usually  randomly 
distributed  over  both  space  and  frequency.  Consequently,  if  the  resolution 
cell  being  tested  contains  only  clutter  (i.e.,  does  not  contain  a  target), 
the  specified  threshold  level  may  still  be  exceeded,  resulting  in  a  false 
target  report,  or  false  alarm.  Clearly,  the  higher  the  setting  of  the 
threshold,  the  lower  the  false  alarm  rate,  but  also  the  lower  will  be  the 
probability  of  detecting  any  given  target.  However,  because  the  number  of 
resolution  cells  to  be  tested  in  a  given  radar  scene  usually  greatly  exceeds 
the  expected  minber  of  targets  in  that  scene,  it  is  important  to  set  the 
threshold  high  enough  to  keep  the  corresponding  false  alarm  rate  at  very  low 
value  -  usually  on  the  order  of  10~3  or  less. 

The  distribution  of  clutter  returns  with  respect  to  spatial  location 
has  been  shown  to  be  log-normal  in  nature  for  many  different  kinds  of  com¬ 
monly  encountered  ground  clutter  [2], [3].  This  distribution  is  a  two-param¬ 
eter  distribution,  and  hence  is  uniquely  characterized  by  both  a  mean  value 


and  a  variance.  The  experimental  results  presented  in  Chapter  5  indicate 
that  the  log-normal  spatial  amplitude  distribution  for  any  specified  type  of 
ground  clutter  (i.e. ,  level  fields,  scrub,  dirt,  deciduous  forest,  conifer 
forest,  etc.)  can  be  uniquely  specified  by  these  two  parameters.  However, 
it  should  be  noted  that,  up  to  the  present  time,  most  quantitative  charac¬ 
terizations  of  clutter  have  been  in  terms  of  only  a  single  parameter,  o0, 
the  mean  value  of  clutter  reflectivity  per  unit  area.  The  mean  value  of  the 
clutter  return  received  by  a  radar  is  proportional  to  Ag^,  where  A  >  tne 
projected  area  of  the  radar's  range-azimuth  resolution  cell  on  the  ci  _cr 
surface.  However  the  value  of  oo  f°r  a  g*ven  clutter  type  generally 
no  indication  of  the  spatial  non-uniformity  of  the  clutter.  This  latter 
quantity  is  characterized  by  the  variance  (or  equivalently  the  standard 
deviation)  of  the  spatial  clutter  amplitude  distribution.  In  this  report 
this  spatial  variance  has  been  given  the  designation  ac2.  As  will  be 
seen,  it  is  usually  convenient  to  work  in  terms  of  the  distribution  of  the 
logarithm  of  clutter  amplitudes.  This  has  the  fortuitous  consequences  that 
1)  the  resultant  distribution  is  normal  (for  log-normally  distributed  clut¬ 
ter)  and  2)  0c  can  be  expressed  in  units  of  dB. 

Therefore,  to  accurately  characterize  the  spatial  amplitude  proba¬ 
bility  distribution  of  a  given  clutter  type,  both  the  mean  reflectivity  0o 
and  the  spatial  standard  deviation,  oc,  oust  be  specified.  Although  many 
measurements  of  go  have  previously  been  reported  in  the  llteratui?  [4], 
(8],  measurements  of  oc  do  not  appear  to  be  as  common.  However,  the  per¬ 
formance  of  an  amplitude  CFAR  detector  in  clutter  depends  upon  both  of  these 
parameters,  [1],  particularly  oc*  The  data  presented  in  Chapter  5 
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indicates  that  for  many  common  types  of  homogenous  ground  clutter  oc 
usually  ranges  between  0.2  and  2.5  dB. 

In  the  context  of  this  report  homogenous  clutter  is  defined  as  clut¬ 
ter  that  is  all  from  the  same  statistical  population  (i.e.,  from  a  region  of 
all  grass,  all  trees,  all  scrub,  etc.),  and  hence  is  characterized  by  a 
unique  set  of  the  parameters:  (o0,  <jc) .  Note  that  this  definition  does 
not  require  that  the  clutter  be  spatially  uniform.  Uniform  clutter  is 
defined  as  clutter  for  which  oc“®»  an<*  *-8  therefore  a  special  case  of 
homogenous  clutter.  Examples  of  uniform  clutter  might  be  an  airport  runway, 
an  empty  parking  lot,  a  level  mown  golf  course,  etc.  It  should  be  noted 
however  that  radar  returns  from  uniform  clutter  will,  in  general,  still  be 
randomly  distributed  with  respect  to  radar  frequency*  Finally  non- 
homogenous  clutter  is  defined  as  mixtures  of  different  kinds  of  homogenous 
clutter  within  a  local  area. 

In  addition  to  being  randomly  distributed,  ground  clutter  is  usually 
also  highly  non-homogenous  over  a  given  radar  surveillance  area.  For 
example,  a  typical  area  may  Include  patches  of  forest,  fields,  roads,  rocky 
hillsides,  etc.,  and  the  statistics  of  the  clutter  distributions  for  each  of 
these  patches  will  be  quite  different  -  i.e.,  the  statistics  (average  value 
and  variance)  of  the  distribution  of  the  returns  from  a  level,  mown  field 
will  generally  be  much  smaller  than  that  for  retwns  from  tree  cover  or  from 
rough  terrain.  Consequently,  if  the  false  alarm  rate  is  to  be  maintained 
constant  over  a  radar  surveillance  area  that  contains  a  variety  of  clutter 
types  such  as  forest  and  fields,  the  CFAR  threshold  must  be  adaptively  set 
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with  respect  to  the  statistics  of  the  local  clutter  in  the  vicinity  of  each 
resolution  cell  being  tested. 


1. 1  Computation  of  CFAR  Detection  Threshold 

Ideally,  if  the  mean  value,  yjj,  and  variance,  ,  Qf  the 
measured  clutter  returns  were  known  exactly,  for  each  local  region  of  the 
radar  surveillance  area,  it  would  be  possible  to  compute  the  threshold  set¬ 
ting  required  to  maintain  a  constant  false  alarm  rate  in  that  region  as 
follows: 


T  “  ux  + 


where  px  and  ox  are  the  mean  and  standard  deviation  of  the  local  clutter 
returns  processed  by  the  CFAR  detector.  The  factor  k  in  the  above  expres¬ 
sion  is  called  the  threshold  multiplier  and  is  uniquely  determined  by  the 
specified  value  of  clutter  false-alarm  rate,  Pfa.  For  typical  false  alarm 
rates  between  lO-^  to  10“^,  the  value  of  k  ranges  between  about  3  and  5. 

The  term  k<jx  is  of  special  significance  since  it  effectively 
represents  the  threshold  offset  relative  to  the  mean  clutter  level.  When 
frequency -averaging  is  used,  yx  and  ox  are  the  mean  and  standard  devia¬ 
tions  of  the  frequency  averages  of  the  clutter  returns.  ox  is  shown  to  be 
a  function  of  both  the  intrinsic  spatial  non-uniformity  of  the  local  clut¬ 
ter,  i.e.,  <jc»  and  the  number  of  frequencies  over  which  the  clutter 
returns  are  averaged.  ox  decreases  at  a  rate  somewhat  less  than  the 
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square-root  of  the  nimber  of  frequencies  over  which  the  returns  are 

averaged,  asymptotically  approaching  the  value  <jc  in  the  limit.  As 
mentioned  previously,  oc  is  equal  to  the  standard  deviation  of  the 
intrinsic  spatial  distribution  of  mean  clutter  amplitudes  (averaged  over 
frequency),  and  is  a  function  of  clutter  type.  Thus  the  value  of  cfc  for 
tree  cover  or  rocky  terrain  will  be  larger  than  that  for  asphalt  pavement  or 
^  a  mown  field. 

For  reasons  to  be  given  later,  it  is  convenient  to  express  the  clut¬ 
ter  statistics  px,  a*  and  oc»  and  the  threshold  level  T  in  terms  of 

dB.  In  terms  of  these  units,  the  value  of  <jc  will  range  between  about  0.5 
and  3  dB  for  the  different  types  of  ground  clutter  encountered  in  a  typical 

radar  scene.  The  corresponding  values  of  ox  will  range  from  about  1.0  dB 

to  7.0  dB,  depending  upon  the  number  of  frequencies  averaged.  As  mentioned 
earlier,  typical  values  of  the  threshold  multiplier,  k,  range  between  3  and 
5,  depending  upon  the  specified  value  of  false  alarm  rate.  Consequently, 
the  resulting  threshold  offset,  kox,  will  usually  range  between  a  low  of 
about  2.0  dB  to  a  high  of  about  30  dB  depending  upon  local  clutter  condi¬ 
tions,  the  number  of  frequencies  averaged,  and  the  specified  value  of  false 
alarm  rate  Pfa.  The  knowledge  of  the  clutter  standard  deviation,  o^,  is 
therefore  especially  critical  in  confuting  the  required  threshold  setting. 

In  practice  however,  values  of  ^  and  ox  of  the  local  clutter 
will  never  be  known  exactly,  but  must  be  estimated  from  the  sample  statis¬ 
tics  of  returns  from  the  local  clutter  in  the  region  being  tested.  As  men¬ 
tioned  above,  because  of  the  different  types  of  clutter  that  can  appear 
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in  a  typical  radar  surveillance  area(  these  local  sample  statistics  may  vary 
widely  over  a  given  radar  scene.  Hence  the  actual  threshold  setting,  T, 
used  for  testing  the  return  from  a  given  cell  is  computed  as  a  function  of 
the  sample  mean,  3,  and  sample  variance,  Vx,  of  the  frequency  averages  of 
returns  from  the  local  clutter  in  the  neighborhood  of  the  cell  being  tested. 

T  «  X  +  k/V^  (1.2) 

Since  the  sample  mean,  1i,  and  sample  standard  deviation,  /Vx,  are  them¬ 
selves  random  variables  -  T  will  also  be  a  random  variable  rather  than  just 
a  deterministic  function  of  the  threshold  multiplier  k.  In  effect  T  and 
Vx  are  merely  estimates  of  the  (unknown)  values  of  the  actual  mean  and 
variance  of  the  local  clutter,  and  <jx^. 

1.2  Effect  of  Frequency  Averaging 

As  can  be  seen  from  equation  1.2,  the  computed  threshold  setting  is 
critically  dependent  upon  the  clutter  sample  variance  Vx.  It  is  shown 
herein  that  if  T  and  Vx  are  computed  from  frequency-averages  of  the  clut¬ 
ter  samples,  rather  than  from  clutter  samples  measured  at  only  a  single  fre¬ 
quency,  the  expected  value  of  5 T  remains  unchanged,  but  that  of  Vx  will 
decrease.  This  is  an  Important  result  since  it  means  that  frequency -averag¬ 
ing  of  the  returns  prior  to  performing  the  CFAR  test  lowers  the  threshold 
setting  required  to  maintain  a  given  value  of  false  alarm  rate  (fixed  k) , 
with  the  result  that  target  detectability  is  improved. 

The  reason  for  this  phenomenon  is  as  follows.  The  amplitude  of 
clutter  returns  varies  randomly  with  respect  to  both  space  and  frequency. 
Specifically,  the  amplitude  of  the  clutter  return  from  any  given  resolution 
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cell  will  be  a  random  function  of  frequency.  The  mean  value,  averaged  over 
frequency,  of  the  returns  from  that  cell  then  represents  the  Intrinsic 
amplitude  of  the  clutter  in  that  cell.  However  these  mean  values  with 
respect  to  frequency  are  in  turn  random  variables  with  respect  to  spatial 
location  (in  range  and  azimuth) ,  and  are  usually  log-normally  distributed 
over  space.  Frequency-averaging  reduces  that  component  of  variability  of 
the  clutter  returns  due  to  frequency.  As  a  result,  the  expected  value  of 
the  sample  variance,  Vx,  that  is  computed  from  frequency-averages  of  clut¬ 
ter  measurements,  decreases  toward  a  lower  bound  equal  to  the  intrinsic 
spatial  variance  of  the  clutter,  <jc^»  88  t^le  nimber  of  frequencies  over 

which  the  clutter  returns  are  averaged  is  increased.  The  more  uniform  the 
clutter,  the  more  effective  frequency-averaging  will  be  in  reducing  Vx. 

The  effect  of  frequency-averaging  on  the  computation  of  the  clutter 
sample  statistics  and  the  resulting  setting  of  the  CFAR  detection  threshold 
is  the  subject  of  this  report.  Of  particular  interest  is  the  effect  on  the 
sample  standard  deviation  since  this  directly  affects  the  required 

level  of  the  CFAR  detection  threshold.  It  is  shown  that  the  expected  value 
of  for  single  frequency  clutter  samples  ranges  from  about  5.5  dB  to 

6.5  dB  for  most  common  types  of  ground  clutter.  As  the  number  of  frequen¬ 
cies  over  which  the  clutter  returns  are  averaged  increases,  the  expected 
value  of  /7^~  approaches  the  value  of  the  standard  deviation  of  the 
clutter' 8  underlying  spatial  distribution  oc  “  8  value  that  is  generally 


between  about  0  and  3  dB ,  depending  upon  clutter  type.  Furthermore,  most  of 
the  attainable  reduction  in  /V^*  occurs  when  averaging  over  between  4  and  8 
frequencies.  Beyond  8  frequencies  a  point  of  diminishing  returns  is  rapidly 
reached  with  regard  to  further  reducing  Averaging  over  4  frequen¬ 
cies,  reduces  the  threshold  offset  (in  dB),  required  to  maintain  a  given 
false  alarm  rate  by  between  35%  and  50%,  depending  upon  clutter  type,  rela¬ 
tive  to  that  required  for  processing  returns  measured  at  a  single  frequen¬ 
cy.  At  8  frequencies,  this  reduction  increases  to  between  40%  and  65%.  The 
results  presented  herein  are  substantiated  by  radar  measurements  of  returns 
from  a  variety  of  different  types  of  ground  clutter  that  were  made  at  Ku 
Band  using  range-azimuth  resolution  cell  sizes  of  36m^  and  108m^. 

It  should  be  pointed  out  that  the  effect  of  frequency-averaging  on 
reducing  the  computed  sample  variance  of  the  clutter,  and  lowering  the 
resulting  CFAR  threshold  level,  is  not  the  whole  story  with  regard  to  the 
effect  of  frequency-averaging  on  improving  fixed  target  detection  perfor¬ 
mance.  Fixed  target  detection  performance  is  characterized  by  the  probabil¬ 
ity  of  detection,  P<j,  realized  at  a  specified  false-alarm  rate,  Pfa,  for 
a  target  whose  return  has  a  given  value  of  target-to-clutter  signal  ratio 
T/C.  Clearly,  as  stated  earlier,  the  lowering  of  the  CFAR  threshold  that 
results  from  frequency-averaging  of  the  clutter  returns  improves  fixed  tar¬ 
get  detectability.  However,  for  a  fixed  target  aspect  angle,  the  target 
returns  are  also  randomly  distributed  as  a  function  of  frequency,  with  a 
standard  deviation  of  about  5.6  dB.  Therefore,  the  standard  deviation  of 
the  frequency-averaged  target  returns  also  decreases  as  the  number  of 


frequencies  over  which  the  returns  are  averaged  is  Increased.  This 
reduction  in  the  standard  deviation  of  the  frequency-averaged  target  return 
provides  an  additional  increase  in  detection  probability,  provided  that  the 
average  value  of  T/C  (for  a  fixed  aspect  angle)  is  greater  than  k<jfc. 
Consequently,  the  improvement  in  fixed  target  detection  performance  due  to 
the  effect  of  frequency-averaging  is  somewhat  greater  than  that  due  to  the 
lowering  of  the  CFAR  threshold  alone,  as  presented  herein.  It  can  be  shown 
[6]  that  the  maximum  number  of  different  frequencies  over  which  returns  can 
be  averaged  before  the  point  of  diminishing  returns  is  reached  with  respect 
to  Improvement  in  fixed  target  detection  performance  is  actually  about  16 
rather  than  the  range  of  4-8  indicated  herein.  However,  a  further  treatment 
of  this  topic  is  beyond  the  scope  of  this  report. 

1.3  Organization  of  Report 

The  report  is  organized  as  follows.  Chapter  2  describes  the  two- 
parameter  CFAR  target  detection  algorithm,  using  frequency-averaged  clutter 
samples  taken  from  a  two-dimensional  array  of  range-azimuth  resolution 
cells.  The  term  "two-parameter"  refers  to  the  fact  that  the  CFAR  threshold 
is  computed  as  a  function  of  both  the  sample  mean  and  sample  variance  of  the 
clutter. 

Chapter  3  describes  the  evaluation  of  false-alarm  probability,  for 
the  two-parameter  CFAR  target  detection  algorithm  described  in  Chapter  2. 
The  analysis  is  based  upon  the  assumption  that  the  frequency-averaged  clut¬ 
ter  samples  processed  by  the  CFAR  algorithm  of  Chapter  2  are  normally 
distributed.  This  assumption  is  a  very  good  approximation  to  reality  when 
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the  clutter  In  the  CFAR  stencil  is  statistically  homogeneous,  with  a 
log-normal  amplitude  distribution  over  space,  and  is  logarithmically 
detected.  Since  the  threshold  is  based  upon  the  computed  sample  statistics, 
the  resultant  false  alarm  probability  is  governed  by  Student's  “t" 
distribution,  as  described  in  Ref.  [1]. 

Chapter  4  describes  the  effect  of  frequency  averaging  on  the  clutter 
sample  statistics  X  and  Vx  that  were  defined  in  Chapter  2.  The  principal 
effect  of  frequency  averaging  is  to  reduce  the  expected  value  of  the  sample 
variance,  Vx.  The  expected  value  of  the  frequency-averaged  clutter  sample 
variance  Vx  is  shown  to  consist  of  two  components  -  the  first  is  the  vari¬ 
ance  of  the  intrinsic  spatial  distribution  of  the  clutter  amplitudes, 

oc  ,  the  second  is  due  to  the  variability  of  the  clutter  returns  from  a 
given  resolution  cell  as  a  function  of  frequency.  It  is  this  second  compo¬ 
nent  of  Vx  that  is  reduced  by  averaging  the  clutter  samples  from  each 
resolution  cell  over  a  number  of  different  frequencies. 

In  Chapter  5,  experimentally  measured  clutter  statistics  for  a  vari¬ 
ety  of  different  types  of  ground  clutter  are  presented.  The  experimentally 
observed  reduction  in  clutter  sample  variance  due  to  frequency-averaging  is 
compared  against  the  reduction  predicted  by  the  analysis  of  Chapter  4.  Very 
good  agreement  between  the  experimentally  measured  and  analytically  predic¬ 
ted  results  is  obtained. 
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2.0  THE  TWO-PARAMETER  CFAR  ALGORITHM 


Consider  the  two-dimensional  array  of  range-azimuth  resolution  cells 
Illustrated  in  Fig.  2.1.  Let  p  be  the  size  of  this  array  in  range  resolu¬ 
tion  cells  and  q  the  azimuthal  extent  of  the  array  in  antenna  beam  posi¬ 
tions.  It  is  assumed  that  the  cells  in  this  array  are  non-overlapping  with 
respect  to  both  range  and  azimuth  so  that  the  returns  from  the  individual 
cells  in  the  array  are  all  statistically  independent.  Normally  p  and  q  are 
odd  numbers  so  that  there  will  be  a  unique  cell  at  the  center  of  this 
array.  All  of  the  cells  in  this  array,  with  the  exception  of  the  center 
cell,  are  assumed  to  contain  local  ground  clutter  returns.  The  size  of  the 
array  is  assumed  to  be  small  enough  so  that  all  the  ground  clutter  within 
the  array  is  homogenous  -  i.e. ,  is  from  the  same  statistical  population. 
Thus  the  effect  of  mixed  types  of  ground  clutter  or  clutter  discontinuities 
such  as  tree  lines  are  not  considered  here.  The  center  cell  is  assumed 
either  to  contain  clutter  from  the  same  population  as  the  surrounding  neigh¬ 
bors,  or  a  target. 

The  presence  of  a  target  in  the  center  cell  of  this  array  can  be 
detected  by  comparing  the  amplitude  of  the  return  from  that  cell  against  a 
threshold  level  that  is  adaptively  computed  as  a  function  of  the  sample  mean 
and  variance  of  the  return  amplitudes  in  remaining  cells  in  the  array, 
excluding  the  cell  under  test.  If  the  return  in  the  center  cell  exceeds 
this  threshold,  a  target  is  declared  to  exist  in  that  cell.  The  amplitudes 
of  the  clutter  cells  in  the  stencil  surrounding  the  cell  under  test  will 
generally  be  randomly  distributed.  The  mean  and  variance  of  their 
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Fig.  2.1  Diagram  of  a  typical  2-dimenslonal  CFAR  stencil,  Illustrating 
the  threshold  detection  process. 


distribution  will  be  a  function  of  the  type  of  clutter  in  those  cells. 


Hence  the  threshold  level  must  be  set  high  enough  so  that  the  probability  of 
a  clutter  return  in  the  center  cell  exceeding  this  threshold  and  being 
declared  a  target  will  be  maintained  at  a  specified  small  value,  usually  on 
the  order  of  10“^  or  less. 

The  amplitudes  of  the  returns  received  from  the  clutter  cells  are 
randomly  distributed  with  respect  to  both  space  and  carrier  frequency  of  the 
radar  signal.  The  effect  of  the  variability  of  the  clutter  amplitude 
returns  as  a  function  of  frequency  may  be  reduced  by  averaging  the  returns 
from  a  given  resolution  cell  over  several  different  frequencies  prior  to 
computing  the  sample  mean  and  variance  of  the  clutter  in  the  stencil  over 
space.  The  advantage  of  using  the  frequency-averaged  returns  from  each 
resolution  cell  in  the  CFAR  test,  rather  than  the  returns  measured  at  only  a 
single  frequency,  is  that  the  threshold  level  required  to  maintain  a  given 
false  alarm  rate  is  lowered,  thereby  increasing  the  probability  of  target 
detection  at  any  given  value  of  false  alarm  rate.  The  probability  of  target 
detection  depends  solely  on  the  target  amplitude  statistics  and  the  thresh¬ 
old  level,  and  is  independent  of  the  clutter  level.  Hence,  the  lower  the 
threshold,  the  higher  the  probability  of  target  detection. 

Let  N“pq-1  be  the  number  of  independent  resolution  cells  in  the  two- 
dimensional  stencil  of  Fig.  2.1,  excluding  the  center  cell.  For  conven¬ 
ience,  let  the  individual  cells  in  this  stencil  be  identified  by  a  one-di¬ 
mensional  index,  i*l,...N.  Let  the  total  number  of  independent  frequencies 


at  which  the  returns  from  each  individual  cell  are  measured  by  J,  and  the 
corresponding  frequency  index  be  j-l,...J. 


The  returns  from  individual  resolution  cells  in  a  given  azimuthal 


beam  cell  can  be  measured  nearly  simultaneously  over  a  number  of  different 
frequencies  by  transmitting  J  pulses  at  a  fixed  antenna  beam  position.  For 
any  reasonable  number  of  frequencies,  the  total  time  required  to  collect 
these  J  pulse  returns  is  usually  only  a  few  milliseconds.  Hence  even  if  the 
radar  is  moving  relative  to  the  ground,  the  motion  of  the  radar  during  this 
interval  is  usually  only  a  very  small  fraction  of  the  resolution  cell  size, 
so  that  for  all  practical  purposes  the  returns  from  all  J  pulses  at  a  given 
beam  position  can  be  assumed  to  be  from  the  same  set  of  physical  resolution 
cells  on  the  ground.  After  each  set  of  J  pulse  returns  have  been  collected 
at  a  given  beam  position,  the  beam  position  is  incremented  by  an  amount 
approximately  equal  to  the  antenna  beam  width  and  the  process  repeated.  The 
data  coiqprislng  the  CFAR  stencil  of  Fig.  2.1  is  then  taken  from  the  qJ  pulse 
returns  that  were  collected  in  covering  q  adjacent  beam  positions.  Let  the 
amplitude  return  from  resolution  cell  i  (i“l , . . .N-pq-1) )  collected  at 
frequency  j  be  designated  as  yjj.  Then  the  frequency-averaged  returns 
from  each  resolution  cell  are  formed  as  follows: 


x 


i 


1 

7 


(2.1) 


The  frequency  average  x^  then  represents  the  best  available  estimate  of 
the  Intrinsic  amplitude  level  of  the  scatterers  in  resolution  cell  i  of  the 
stencil.  Furthermore,  designate  i*0  as  the  index  of  the  center  cell  of  the 


stencil.  Assume  that  the  N  cells  in  the  stencil,  outside  of  the  central 


cell,  all  contain  local  clutter  of  a  given  type.  The  sample  mean  and 
variance  of  the  cell  frequency  averages,  x*,  over  the  stencil  may  be  then 
used  as  estimates  of  the  mean  and  variance  of  the  spatial  probability 
distribution  of  the  local  clutter.  These  sample  statistics  are  computed  as 
follows.  Clutter  Sample  Mean: 


Clutter  Sample  Variance: 


V 

x 


(2.3) 


It  should  be  noted  that  the  value  of  the  sample  mean,  X,  is  the  same  regard¬ 
less  of  whether  this  mean  Is  taken  of  the  frequency  averages  x*,  or  of  the 


raw  frequency  returns  yjj  over  both  frequency  and  space.  However,  this  Is 


not  true  of  the  sample  variance  Vx. 

l  N  N 

vx  +  ir  £  £  (yA1 

J  i*l  j-l  1J 


For  example: 

-X)2 


Note  that  the  return  xQ,  from  the  center  cell  of  the  stencil  Is  excluded 
from  the  above  computation.  To  complete  the  description  of  the  CFAR  detec¬ 
tion  process,  a  threshold  value  Is  computed  Individually  for  each  clutter 
stencil  as  follows: 

T  -T  +  (2.4) 
x 


Where  is  the  square-root  of  equation  (2.3).  The  frequency  average  of 
the  center  cell  x0  is  then  compared  against  this  threshold: 


>  T? 


x 

o 

If  x0  is  greater  than  T,  a  target  is  assumed  to  exist  in  the  center  cell. 
However,  even  if  the  center  cell  contain*,  only  clutter  -  i.e.,  if  xQ  is 
from  the  same  statistical  population  as  its  neighbors  x^,  i“0,  there  is 
still  a  finite  probablity  that  xQ>T,  resulting  in  a  false  alarm.  The 
probability  of  such  a  false  alarm  occurring  is  a  function  of  the  value  of  T 
and  the  probability  distribution  of  the  x^,  as  discussed  in  detail  in  the 
next  section. 
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3.0  EVALUATION  OF  FALSE  ALARM  PROBABILITY 

The  variable  used  in  the  CFAR  test  is  the  frequency  average  of  the 
radar  returns,  x^.  As  is  shown  in  Appendix  B,  if  the  x^  represent  fre¬ 
quency  averages  of  log-detector  outputs  in  response  to  radar  returns  from 
log-normally  distributed  clutter,  the  distribution  of  the  x^  will  be 
approximately  normal,  with  a  mean  and  variance  that  are  a  function  of  the 
clutter  statistics.  This  situation  is  illustrated  in  Fig.  3.1. 

Then  for  any  given  threshold  setting,  T,  the  false  alarm  probability 
Pfa  can  be  computed  as  the  probability  of  x  >T  and  is  represented  by  the 
area  under  the  (approximately  normal)  probability  density  function  f(x)  of 
the  frequency-averaged  clutter  amplitudes  between  T  and  infinity  as 
illustrated  in  Fig.  3.1. 

OO 

Pfa(T)  "  T /  f(x)dx)  (3.0 

Furthermore,  if  the  density  function  f(x)  is  normal  and  the  thresold  T  is 
computed  as  an  offset  from  the  mean,  yx,  of  f(x)  as: 

T  -  Ux  +  k0x  (3,2) 

Then  it  can  be  easily  shown  that  Pfa,  as  defined  by  equation  3.1,  is 
merely  a  function  of  the  threshold  multiplier  k,  i.e., 

P,  (k)  -  f“f(x)dx  -  Erfc(k)  (3.3) 

u  +k(j 

MX  X 

In  equation  3.2  above,  the  mean  and  variance  of  x,  i.e.,  yx  and 

2 

Ox  are  assumed  to  be  known  parameters  -  i.e.  ,  constants.  In  reality, 
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Fig.  3.1  Probability  distributions  of  target  and  clutter  returns  for  a 
given  value  of  target-to-clutter  ratio,  indicating  effect  of  detection 
threshold  setting  on  target  detection  and  clutter  false  alarm  probabili¬ 
ties.  Effect  of  frequency-averaging  on  narrowing  of  the  distributions  is 
also  illustrated. 


however,  these  parameters  are  unknown,  and  T  can  only  be  computed  from 
estimates  of  px  and  ox  ~  i»e.,  from  the  sample  statistics  as  in  equation 
2.4.  As  a  result,  T  Itself  becomes  a  random  variable.  It  can  be  shown  [1] 
that  when  T  is  computed  from  the  sample  statistics  X  and  Vx,  (equation 
2.4)  and  the  are  normally  distributed,  the  probability  of  false  alarm 

is  given  by  the  following  integral: 


Pfa(N,k)  -  J  SN_1(Z)dZ  (3.4) 


where  Sj}_i  is  Student's  distribution  for  N-l  degrees  of  freedom,  and  N  is, 
as  before,  the  number  of  resolution  cells  in  the  stencil  used  to  compute  the 
sample  mean  and  variance,  X  and  Vx.  Specifically,  Student's  distribution 
Sn_i(Z)  is  the  distribution  of  the  random  variable: 


Z  - 


where  X  and  Vx  are  as  defined  by  equations  2.2,  2.3  and  the  sample  x0  is 
from  the  same  population  as  the  x^,  but  is  not  included  in  the  computation 
of  X  and  Vx.  As  might  be  expected,  as  N+  °%  Student's  distribution 
approaches  a  zero-mean  normal  distribution  of  unit  variance.  Indeed  for  all 
practical  purposes,  when  N  >40  the  value  of  Pfa  computed  from  equation  3.4 
approaches  that  of  equation  3.3,  for  T  defined  by  equation  2.4. 

In  summary,  when  the  x^  are  normally  distributed,  and  T  is  com¬ 
puted  as  a  function  of  the  sample  mean  and  variance  in  accordance  with  equa¬ 
tion  2.4,  the  false  alarm  probability  will  be  a  unique  function  of  the 
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threshold  multiplier,  k,  and  the  nunber  of  resolution  cells,  N,  over  which 
the  sample  statistics  are  computed.  However,  in  the  sequel,  the  effect  of  N 
on  Pfa  will  not  be  considered  -  i.e.,  N  will  be  assumed  to  be  fixed. 

4.0  EFFECT  OF  FREQUENCY -AVERAGING  ON  THE  CFAR  THRESHOLD  SETTING 

The  consequence  of  the  previous  argument  that  Pfa  is  only  a  func¬ 
tion  of  the  threshold  multiplier,  k,  (the  effect  of  N  being  neglected  since 
N  is  usually  fixed)  is  that  in  order  to  maintain  a  specified  Pfa,  the 
absolute  value  of  the  threshold  setting  T  will  be  proportional  to  both  the 
sample  mean,  X,  and  the  square  root  of  the  sample  variance,  /Vx,  in  accor¬ 
dance  with  equation  2.4.  However  in  order  to  maximize  the  probability  of 
target  detection  P<j,  the  threshold  should  be  kept  as  low  as  possible,  con¬ 
sistent  with  maintaining  the  specified  value  of  Pfa.  The  question  at  hand 
then  is  how  can  frequency  averaging  of  the  data  be  used  to  lower  the  value 
of  the  threshold  setting,  T,  without  raising  Pfa? 

The  major  benefit  provided  by  frequency  averaging  is  that,  for  most 
types  of  clutter,  frequency  averaging  significantly  reduces  the  expected 
value  of  the  clutter  sample  variance,  Vx.  This  allows  the  CFAR  threshold 
level  required  to  maintain  a  given  false  alarm  rate  to  be  decreased  corre¬ 
spondingly,  as  per  equation  2.4,  with  the  result  that  higher  detection  prob¬ 
abilities  can  be  realized  for  a  given  value  of  clutter  false  alarm  rate. 


It  Is  shown  in  Appendix  A  that  the  expected  value  of  the  sample 


variance,  Vx,  of  the  cell  frequency  averages,  X^,  is  given  by  the 
following  expression. 

EIV  ■  T  (7  siv.J  +  °c2)  «•*> 

In  the  above  expression  is  the  variance  of  the  spatial  distribution 

2 

of  the  intrinsic  clutter  amplitudes,  and  E  [oy/m  ]  is  the  expected  value  of 

the  variance  of  the  clutter  amplitude  returns  received  from  a  resolution 
cell  containing  clutter  of  intrinsic  amplitude  level  mx,  and  is  a  measure 
of  the  variability  of  the  returns  from  that  cell  over  frequency.  J  is  th'j 
number  of  frequencies  averaged,  and  N  is  the  number  of  independent  resolu¬ 
tion  cells  in  the  stencil.  N  should  be  kept  small  enough  so  that  the  sten¬ 
cil  contains  only  homogenous  clutter  of  variance  <jc^»  but  large  enough 
so  that  the  sample  statistics,  X  and  Vx,  are  good  estimates  of  the  true 
mean  and  variance  of  the  local  clutter.  A  reasonable  compromise  is  to 

choose  a  stencil  having  between  20  and  40  resolution  cells.  The  value  of 
2 

oc  is  difficult  to  measure  and  will  never  be  known  exactly.  As  previ¬ 
ously  discussed,  it  is  a  characteristic  of  the  local  clutter,  and  may  vary 

greatly  within  any  given  radar  scene.  However,  as  can  be  seen  from  equation 
4.1,  oc^  can  be  estimated  from  Vx  when  J  is  large  -  i.e.  ,  when  the 
clutter  returns  from  a  given  ground  patch  are  averaged  over  a  large  number 


of  frequencies 


Effect  of  Logarithmic  Detection 


If  the  clutter  returns  are  logarithmically  detected,  the  output  of 

the  log  detector,  being  proportional  to  the  logarithm  of  the  signal  voltage, 

will  be  linear  in  terms  of  dB  relative  to  some  arbitrary  reference.  Hence, 

the  averages,  x^,  of  the  log  detector  outputs  can  also  be  expressed  in  dB 

relative  to  this  reference.  Clearly,  the  choice  of  the  reference  level  only 

effects  the  origin  of  the  dB  scale.  Hence,  although  the  mean  value,  T,  of 

the  x±  will  be  dependent  upon  this  choice  of  origin,  the  variance,  Vx, 

2  2 

of  the  x<  will  not.  Therefore  the  terms,  a  and  E[<j  ,  ] ,  in  equation 

x  c  y/ mx 

4.1  above  can  be  expressed  in  dB^,  and  are  independent  of  the  choice  of 

origin  for  the  dB  scale. 

Furthermore,  it  is  shown  in  Appendix  C  that  for  the  case  of  a 


logarithmic  detector,  the  quantity  o 


31  dB  is  a  constant, 


independent  of  mx,  whose  value  is: 


31  dB. 


4.2  Clutter  Sample  Variance  as  a  Function  of  Frequencies  Averaged 
and  Intrinsic  Clutter  Variance 

Substituting  this  into  equation  4.1  gives  the  following  result  for  the 
expected  value  of  the  stencil  variance: 


“V-W  +°c2< 


(4.2) 


The  effect  of  frequency  averaging  on  variance  reduction  is  now  clear.  Fre¬ 
quency  averaging  only  reduces  that  component  of  the  total  observed  clutter 


variance  -  i.e.,  31  dB^  -  that  is  attributable  to  variation  of  the 

amplitude  of  the  clutter  returns  with  frequency.  Hence  the  effectiveness  of 
frequency  averaging  depends  upon  the  magnitude  of  the  intrinsic  clutter 
variance,  oc^,  relative  to  31  dB^,  or  equivalently  the  standard  devia¬ 
tion  of  the  clutter  relative  to  5.6  dB.  Clearly  the  smaller  the  intrinsic 
standard  deviation  of  the  clutter,  ac»  the  more  effective  frequency  aver¬ 
aging  will  be  in  reducing  the  expected  value  of  the  stencil  variance,  Vx. 
This  phenomenon  will  be  explored  in  more  detail  and  compared  with  a  variety 
of  experimental  clutter  measurements  in  the  next  section. 

4.3  Effect  of  Frequency-Averaging  on  Threshold  Setting 

As  mentioned  earlier,  the  effect  of  frequency  averaging  on  OAR 
detection  performance  is  to  lower  the  threshold  setting  required  to  maintain 
a  given  false  alarm  rate,  thereby  improving  detection  performance  by 
increasing  the  probability  of  target  detection  without  compromising  false 
alarm  performance.  This  can  be  seen  by  substituting  equation  4.2  for  the 
expected  value  of  Vx  into  equation  2.4  for  the  threshold  setting,  and  is 
also  illustrated  in  Fig.  4.1.  Hence  an  expression  for  the  effect  of 
frequency  averaging  on  the  expected  value  of  the  threshold  setting  required 
to  maintain  a  given  value  of  false  alarm  rate  is  given  by: 

E[T]  -  ux  +  k/^L  (^1  +  oc2J  dB  (4.3) 

Or  equivalently,  the  expected  value  of  the  threshold  setting  relative  to  the 
mean  value  of  the  clutter  is: 
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Fig.  4.1  Probability  distribution  of  log-detected  clutter  amplitudes  with 
and  without  frequency-averaging  illustrating  the  effect  of  frequency¬ 
averaging  on  narrowing  the  distribution  and  reducing  the  threshold  level 
required  to  maintain  a  specified  Pfa* 


(4.3a) 


As  previously  discussed,  the  false  alarm  rate  corresponding  to  this 
threshold  setting  is  determined  by  the  value  of  the  threshold  multiplier  k, 
for  fixed  N.  Note  that  equation  (4.3a)  for  the  relative  threshold  setting 
is  independent  of  the  dB  scale  origin  chosen  for  the  log  detector  outputs. 
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5.0 


EXPERIMENTAL  VERIFICATION 


In  the  preceding  section,  a  quantitative  relation,  equation  (4.2), 
was  obtained  for  the  expected  value  of  the  clutter  sample  variance  as  a 
function  of  the  number  of  frequencies  and  number  of  resolution  cells  aver¬ 
aged,  and  the  variance  of  the  underlying  spatial  amplitude  distribution  of 
the  clutter  in  the  stencil.  Consequently,  equation  4.2  may  be  used  to  pre¬ 
dict  the  effect  of  frequency  averaging  on  the  sample  variance  computed  for 
different  types  of  clutter. 

For  this  purpose,  each  clutter  type  may  be  uniquely  characterized  by 
the  variance,  oc2.  of  its  intrinsic  spatial  amplitude  distribution.  For 
example,  smooth  clutter  such  as  mown,  level  fields,  or  pavement  is  charac¬ 
terized  by  a  spatial  amplitude  distribution  having  a  very  low  value  of 
< jc whereas  tree  cover  or  rough  terrain  would  be  characterized  by  a 
spatial  amplitude  distribution  having  a  much  higher  value  of  0c 2. 
Ideally,  the  values  of  oc2  for  each  clutter  type  can  be  obtained  as  the 
limit  of  the  sample  variance,  over  space,  of  clutter  returns  that  have  been 
frequency-averaged  over  an  infinite  number  of  independent  frequencies. 
Clutter  returns  obtained  at  two  different  frequencies  may  be  shown  to  be 
statistically  independent  if  the  difference  between  the  two  frequencies  is 
greater  than  the  radar  ban<fcfidth  17).  Obtaining  measurements  at  an  infinite 
number  of  frequencies  is  clearly  impractical;  however  in  most  cases  of  prac¬ 
tical  interest,  the  limiting  variance  of  the  spatial  clutter  amplitude  dis¬ 
tribution  will  be  approached  asymptotically,  within  the  margin  of  experimen¬ 
tal  error,  when  clutter  returns  are  averaged  over  more  than  32  frequencies. 
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5.1  The  Variance  Reduction  Factor 


Of  particular  Interest  la  the  reduction  In  the  expected  value  of 
clutter  sample  variance  .  Vx,  that  results  from  averaging  the  clutter 

returns  over  a  specified  number  of  frequencies.  It  Is  the  ratio  of  the 

expected  value  of  Vx,  when  the  clutter  samples  are  averaged  over  J  fre¬ 

quencies,  to  the  expected  value  of  Vx  computed  from  clutter  samples  meas¬ 
ured  at  only  a  single  frequency  (i.e.,  no  frequency  averaging).  This  vari¬ 
ance  reduction  ratio  is  defined  as  p(J,oc)  and  is  obtained  from  equation 

4.2  as  follows: 

E [V  (J)]  31  +  J02  .  2 

p(J,0  )  -  - - -  - - £  where  °c  and  31  are  (5.ls 

C  E[Vx(l)]  J(31+ac2)  units  of  dB2. 

The  functional  dependence  of  this  variance  reduction  factor  on  both  the  num¬ 
ber  of  frequencies,  J,  averaged,  and  upon  the  standard  deviation,  0c,  of 
the  Intrinsic  spatial  clutter  amplitude  distribution  is  specifically  indi¬ 
cated. 

5.2  Description  of  Experimental  Clutter  Measurement  Data 

Experimentally  derived  clutter  statistics  for  a  wide  variety  of  com¬ 
monly  encountered  types  of  ground  clutter  have  been  measured  at  Ku  band  and 
reported  in  [2].  The  results  of  these  measurements  are  summarized  in  Tables 
5.1  and  5.2.  Table  5.1  is  a  summary  of  clutter  measurements  made  at 
Stockbridge,  NY  in  the  summer  of  1977,  using  a  100m2  range-azimuth  resolu¬ 
tion  cell.  The  measurements  were  averaged  over  32  different  frequencies 
spanning  a  bandwidth  of  500  MHz.  Table  5.2  summarizes  the  results  of  a 


TABLE  5,1.  Summary  of  Clutter  Parameters* 
Stockbridge,  NY  -  Late  Summer.  1977 


*Resolution  cell  (9  m,  range)  x  (12  m,  cross 
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similar  sec  of  clutter  measurements  made  at  Fort  Devens,  HA  in  the  summer  of 
1978  using  a  40m2  range-azimuth  resolution  cell  and  16  different  frequen¬ 
cies  spanning  the  same  500  MHz  bandwidth. 

The  smaller  range-resolution  used  at  Fort  Devens  corresponds  to  a 
wider  radar  signal  bandwidth.  Therefore,  in  order  to  maintain  statistical 
Independence  between  the  measurements  made  at  different  frequencies  at  the 
wider  signal  bandwidth,  the  differential  between  adjacent  frequencies  had  to 
be  increased.  Since  the  overall  frequency  agility  bandwidth  of  the  radar 
was  limited  to  500  MHz,  the  required  increase  in  frequency  differential  was 
realized  by  reducing  the  number  of  frequencies  spanning  that  bandwidth  from 
32  to  16. 


The  clutter  statistics  are  represented  in  terms  of  the  clutter 
reflectivity,  y,  which  is  a  dimensionless  quantity  corresponding  to  apparent 
cross-section  of  the  clutter  return  in  m^,  divided  by  the  projection  of 
the  ground  area  covered  by  the  radar's  range-azimuth  resolution  cell  onto  a 
plane  normal  to  the  radar  line-of-sight.  Another  dimensionless  measure  of 
clutter  reflectivity  that  is  often  used  is  0O,  defined  as: 

apparent  radar  cross-section 


o  (.ground  area  ot  radar  resolution  ce 


The  two  parameters  are  related  as  follows: 

Y  “  0o  C8C  0d 

where  e<j  is  the  depression  angle  ol  the  radar-line-of-sight  with  respect 
to  the  ground  plane.  Because  of  this  additional  normalization,  y  is  an 
intrinsic  clutter  parameter  that  is  more  nearly  independent  of  depression 
angle  over  a  wide  range  of  values  of  qh  [4]. 
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As  pointed  out  in  the  last  chapter,  since  the  clutter  amplitude 
returns  are  logarithmically  detected,  they  may  be  represented  directly  in 
terms  of  dB .  Thus  the  above  normalization  only  results  in  a  shift  of  the  dB 
scale  origin  and  has  no  effect  on  the  standard  deviation  (or  variance)  of 
the  clutter  returns,  when  expressed  in  dB.  Therefore  the  previously  derived 
expressions  for  clutter  sample  variance  in  dB  are  still  applicable  when  the 
clutter  returns  are  expressed  in  terms  of  y  (dB). 

The  meaning  of  the  various  columns  of  Tables  5.1  and  5.2  are  defined 
in  Table  1.  However  a  few  amplifying  comments  may  be  helpful.  Column  1  is 
the  mean  value  of  the  log-detected  clutter  returns  in  dB  averaged  over  both 
space  and  frequency.  It  is  computed  using  the  sacond  expression  for  X  in 
equation  2.2.  Column  3  is  the  spatial  mean,  over  the  stencil,  of  the  fre¬ 
quency  averages  of  the  log-detected  clutter  returns  for  each  cell  of  the 
stencil.  It  is  computed  using  the  definition  of  the  frequency  average  given 
by  equation  2.1  and  the  first  expression  in  equation  2.2.  As  would  be 
expected  from  the  equality  of  the  two  expressions  in  equation  2.2,  the 
values  in  these  two  columns  are  identical.  Column  2  represents  the  sample 
standard  deviation  of  the  clutter  returns  measured  at  only  a  single  frequen¬ 
cy,  in  effect  .  Column  4  represents  the  sample  standard  deviation 

of  the  frequency-averages  as  defined  by  equation  2.3,  in  effect  /V^(J). 
Column  5  again  represents  the  mean  value  of  the  clutter,  but  where  the  aver¬ 
aging  is  done  prior  to  logging.  Hence,  the  difference  between  columns  (1) 
and  (3)  on  the  one  hand  and  column  (5)  on  the  other,  is  that  the  values  in 
column  (1)  and  (3)  are  computed  as  the  averages  of  the  log-detected  clutter 


signal  voltages,  whereas  column  (5)  is  computed  as  the  logarithm  of  the 
average  clutter  signal  voltage.  It  can  be  shown  that  the  two  quantities  are 
related  by  the  expression  given  in  footnote  (5)  of  Table  5.1.  In  that 
expression  the  numbers  in  parentheses,  including  the  (2),  refer  to  the 
column  numbers  of  the  table. 

Finally,  coliann  (6)  of  the  tables  is  the  reduction  in  the  sample 
variance  of  the  frequency  averaged  clutter  returns  that  results  from  fre¬ 
quency  averaging,  and  is  computed  as  the  square  of  the  ratio  of  column  (2) 
(the  sample  standard  deviation  for  returns  measured  at  a  single  frequency) 
to  column  (4)  (the  sample  standard  deviation  for  returns  frequency-averaged 
over  J  frequencies,  where  J-32  for  Table  5.1  and  J-16  for  Table  5.2).  Thus 
the  data  in  column  (6)  of  the  tables  corresponds  to  the  reciprocal  of  the 
variance  reduction  factor  p(J,<jc)  defined  by  equation  (5.1). 

5.3  Estimation  of  Intrinsic  Spatial  Clutter  Standard  Deviation 
From  Measured  Clutter  Sample  Statistics 

In  order  to  compare  the  experimentally  measured  variance  reduction 
factors  against  the  predicted  values  from  equation  5.1,  it  is  first  neces¬ 
sary  to  determine  the  value  of  the  standard  deviation  of  the  intrinsic  spat¬ 
ial  distribution  of  clutter  amplitudes,  oc>  (corresponding  to  averaging 
over  an  infinite  number  of  frequencies)  for  each  of  the  clutter  types  listed 
in  Tables  5.1  and  5.2. 

This  process  is  a  bit  tricky  because  only  two  frequency-averaged 
data  points  are  available  for  each  type  of  clutter,  i.e. ,  single-frequency 
(no  averaging)  and  the  averages  taken  over  either  16  or  32  frequencies,  and 
also  because  of  probable  random  measurement  and  calibration  errors  in  the 
data,  as  will  be  discussed  further  below.  Ideally,  an  estimate  of  oc^ 
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would  be  obtained  for  each  clutter  type  by  computing  the  experimentally 

derived  variance  reduction  factor  VX(J)/VX(1)  from  the  data  measured  at 
k  different  frequencies,  at  a  number  m  (>3)  different  values  of  J,  Jin<k> 
so  as  to  obtain  a  sequence  of  values  of  Vx(Jm)/Vx(l )  for  that  type  of 
clutter,  for  m  different  numbers  of  frequencies  averaged.  A  least  squares 
fit  of  equation  5.1  would  then  be  made  to  the  resulting  sequence  with 
oc ^  as  a  parameter.  The  value  of  ac^  that  produced  the  best  fit  of 
equation  (5.1)  to  this  sequence  of  m  values  of  Vx(Jm) /Vx(l )  would  be 
the  best  estimate  of  ac2  for  that  type  of  clutter.  Unfortunately,  the 

values  of  VX(J)  cited  herein  are  only  available  at  two  values  of  J  for 
each  clutter  type:  either  J=16  or  32,  and  trivially  for  J=l.  Although  a  fit 
of  equation  (5.1)  could  be  obtained  that  would  pass  through  the  resulting 
two  data  points,  i.e.  ,  Vx(  1 ) /Vx(  1 )  =1  and  VX(J)/VX(1),  for  each  clut¬ 
ter  type,  the  problem  is  further  complicated  by  the  fact  that  there  is  also 

an  apparent  calibration  error  in  the  data,  as  discussed  below. 

From  equation  (4.2)  it  follows  that: 


N 

TFT 


(5.2) 


The  number  of  resolution  cells  used  to  obtain  the  values  of  VX(J)  tabula¬ 
ted  in  Tables  5.1  and  5.2  ranged  from  200-800.  Hence  for  all  practical  pur¬ 
poses  the  factor  (N/N-l)  »  1.  However  a  calibration  discrepancy  is  indica¬ 
ted  by  the  fact  that  with  one  exception,  the  values  of  /Vx(.  1 )  listed  in 
column  (2)  of  Tables  5.1  and  5.2  are  all  less  than  /TT,  so  that  if  the 
corresponding  values  of  Vx(l)  were  substituted  into 
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equation  5.2,  negative  values  of  oc2  would  be  obtained  -  a  clearly 
unrealistic  result.  Consequently  the  values  of  Vx(l)  indicated  in  Tables 
5.1  and  5.2  appear  to  be  somewhat  low.  The  problem  appears  to  be  due  to 
errors  in  the  calibration  of  the  logarithmic  detector.  Although  the  exact 
cause  of  this  discrepancy  is  unknown,  a  likely  cause  appears  to  be  the  band¬ 
width  limit  of  the  logarithmic  detector.  The  log-detector  is  effectively 
calibrated  at  "D.C."  -  i.e. ,  the  detector  outputs  are  measured  at  constant 
input  signal  levels  that  are  incremented  in  a  step-by-step  manner  in  order 
to  obtain  the  detector  calibration  curve.  However,  if  the  amplitude  varia¬ 
tion  of  the  clutter  returns  between  adjacent  range  resolution  cells,  divided 
by  the  range  resolution  in  units  of  time,  is  greater  than  the  detector  band¬ 
width,  the  magnitude  of  the  cell-cell  fluctuations  in  clutter  amplitude  will 
be  attenuated.  Consequently,  the  variance  of  the  log-detected  returns  will 
be  somewhat  less  than  their  true  variance,  giving  rise  to  the  observed 
behavior.  As  a  result,  in  order  to  use  equation  (5.2)  to  estimate  oc2 
for  the  data  in  Tables  5.1  and  5.2,  the  constant:  31  dB^  will  be  replaced 
by  the  quantity  Vx(l)-0c2,  where  Vx(l)  is  obtained  from  column  2  of 
Tables  5.1  and  5.2.  This  will  force  equality  of  equation  (5.2)  when  J=l, 
regardless  of  the  measured  value  of  Vx(l).  Note  that  the  value  of  Vx(l) 
is  only  weakly  influenced  by  the  value  of  the  clutter  variance  ac2* 
those  cases  where  ©c  <  3  dB .  Hence  assuming  that  N/N-l  »  1,  a  modified 
form  of  equation  (5.2),  adjusted  for  the  effects  of  the  apparent  calibration 


error  is: 


(5.3) 

The  nunber  of  frequencies,  J,  is  either  32  or  16,  depending  upon 
whether  the  data  from  Table  5.1  or  5.2  is  being  adjusted.  The  expected 
clutter  variance  at  the  indicated  mrnber  of  frequencies,  E[VX(J)],  is 
obtained  as  the  square  of  the  standard  deviation  of  the  frequency  averages 
given  in  column  (4)  of  the  tables.  Estimated  values  of  <jc  for  each  clut¬ 
ter  type  in  Tables  5.1  and  5.2  are  obtained  by  substituting  the  observed 
values  of  Vx(l)  and  VX(J)  into  equation  (5.3).  The  results  are  tabulat¬ 
ed  in  Table  5.3  together  with  the  number  of  frequencies  and  the  experimen¬ 
tally  measured  variance  reduction  factor  (the  reciprocal  of  column  (6)  in 
Tables  5.1  and  5.2).  Having  estimated  the  value  of  the  intrinsic  clutter 
variance,  0c,  for  each  clutter  type  listed  in  Tables  1  and  2,  one  can  now 
compare  the  corresponding  experimentally  derived  variance  reduction  factors, 
pU»^c)»  against  the  values  predicted  using  equation  5.1. 

5.4  Comparison  of  Experimental  Data  with  Theoretical  Predictions 

It  is  somewhat  more  useful  to  portray  the  effect  of  frequency¬ 
averaging  in  terms  of  the  reduction  in  standard  deviation  of  the  clutter 
rather  than  the  reduction  in  variance,  since  it  is  the  estimated  standard 
deviation  that  affects  the  detection  threshold  level.  Therefore  the 


A  2 


vx(j)  - 
A 


V  (l)-'a 


_A  2 


/\  9 

Or  solving  the  oc  gives : 


A  2 

°c 


J  Vx(1) 

A  VJ>  -  4t- 


35 


frequency-average  standard  deviation  reduction  factor  /p(J ,  aci »  based  upon 
equation  5.1,  is  plotted  in  Fig.  5.1  as  a  function  of  the  number  of  frequen¬ 
cies  averaged,  J,  for  clutter  standard  deviations,  oc»  ranging  from  0  to 
3.0  dB,  in  steps  of  1.0  dB.  Note  that  the  curve  for  <jc  ■  0  dB  represents 
the  ideal  reduction  of  /p(J,0)  ■  /1/T. 

Also  plotted  in  Fig.  5.1,  for  comparison,  are  the  experimentally 
derived  standard  deviation  reduction  factors  based  upon  the  measurement  data 
tabulated  in  Tables  5.1  and  5.2.  The  estimated  value  of  for  each 
experimental  data  point,  that  was  obtained  using  equation  (5.3)  and  is 
tabulated  in  Table  5.3,  is  also  indicated  in  Fig.  5.1.  The  agreement 
between  the  experimentally  derived  data  points  and  the  theoretical  reduction 
curves  predicted  using  equation  (5.1)  appears  to  be  quite  good  for  all  cases 
considered,  i.e. ,  the  values  of  oc  estimated  for  each  experimentally 
derived  data  point  are  within  less  than  0.5  dB  of  the  theoretical  value  of 
oc  that  would  yield  the  observed  value  of  reduction  at  the  indicated  num¬ 
ber  of  frequencies. 

As  mentioned  previously,  the  experimentally  derived  data  tabulated 
in  Tables  5.1  and  5.2  appear  to  be  subject  to  a  number  of  measurement  and 
calibration  errors.  Such  errors  can  effect  the  estimated  values  of  oc* 
tabulated  in  Table  5.3,  by  as  much  as  0.5  dB.  The  effect  of  these  errors  on 
p(J,oc)  is  believed  to  be  less  because  this  parameter  is  computed  as  the 
ratio  of  two  measurements,  and  systematic  errors  in  the  measurements  would 
tend  to  cancel  in  this  ratio.  Two  sources  of  such  errors  are  known. 
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space)  for  different  clutter  types* 


First,  as  already  mentioned,  there  appears  to  be  a  discrepancy  of 
about  0.5  dB  in  the  single-frequency  clutter  standard  deviations  tabulated 
in  columns  (2)  of  Tables  1  and  2.  Although  this  discrepancy  is  small  and 
appears  to  be  on  the  order  of  the  expected  experimental  error,  it  appears  to 
be  too  systematic  to  be  random.  The  exact  cause  is  unknown,  but  as 
mentioned  previously,  it  is  believed  to  be  due  to  dynamic  miscalibratlon  of 
the  logarithmic  detector  at  high  signal  bandwidths. 

Second,  there  is  an  effective  quantization  error  of  0.1  dB  in  the 
data  of  Tables  1  and  2,  in  that  the  data  were  only  tabulated  with  that 
precision.  However  an  error  of  0.1  dB  in  the  value  of  /EIVX]  listed  in 
column  (4)  of  the  tables  can  result  in  an  error  in  the  estimated  value  of 
nr  obtained  from  equation  5.2  of  as  much  as  0.3  dB,  at  the  lower  values  of 
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APPENDIX  A 


Derivation  of  CFAR  Stencil  Variance  as  a  Function  of  Stencil  Size, 
Number  of  Frequencies  Averaged,  and  Clutter  Statistics 

Consider  the  two-dimensional  range-azimuth  CFAR  stencil  shown  in 
Fig.  A-l,  consisting  of  N*pxq  independent  resolution  cells,  where  p  is  the 
extent  of  the  stencil  in  range  cells,  and  q  is  the  azimuthal  extent  in  beam 
positions.  Let  each  of  the  N  resolution  cells  in  this  stencil  be  identified 
by  an  index  i,  i*l,  ...  N.  Let  returns  be  collected  from  each  cell  in  the 
stencil  at  J  different  frequencies,  sufficiently  separated  so  that  the  cor¬ 
responding  J  returns  from  the  given  cell  are  all  statistically  independent. 
Let  the  amplitude  of  the  return  from  cell  i  for  frequency  j  be  designated 
y^j.  The  distribution  of  the  y^j  is  two-dimensional.  One  dimension 
represents  frequency,  the  other  the  spatial  variation  (i.e.  ,  cell  location) 
in  the  stencil.  The  mean  value  of  the  return  from  each  cell,  over  fre¬ 
quency,  represents  the  intrinsic  amplitude  of  the  clutter  in  that  cell,  and 
is  a  random  variable  with  respect  to  space.  Let  the  mean  value  of  the  fre¬ 
quency  returns  from  cell  i  be  Then  the  two-dimensional  probability 

distribution  function  of  the  returns  can  be  represented  by: 

p(y,z)  »  py/2(y/z)  pz(z)  (A-l) 

Note  that  y  and  z  are  not  independent  rv.;  in  effect  z  is  the  mean  value  of 
the  distribution  Py/z(y/z). 
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Lot  Xj  be  the  frequency-averaged  amplitude  of  the  radar  return  in  the  ith  cell  of  the 
stencil;  let  Xc  be  the  frequency-averaged  amplitude  of  the  return  in  the  center  cell 


Define; 
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Xi 
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Detection  Threshold  T  =  X+k 


Then  if  X^T,  a  target  is  declared  to  be  present  in  the  center  cell 


Fig.  A-l.  Diagram  of  a  typical  2-dimensional  CFAR  stencil,  illustrating 
the  threshold  detection  process. 
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(A. 2) 
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However  as  argued  by  Stovall^],  Py/z(y/z)  and  Pz(z)  are  two  distinct 
identifiable  distributions,  Pz(z)  is  the  spatial  distribution  of  the  mean 
clutter  amplitude  and  is  generally  log-normal,  whereas  Py/z(y/z)  is  the 
amplitude  distribution  of  the  frequency  returns  from  a  given  cell  and  is 
generally  exponential  or  Rayleigh  depending  upon  whether  the  amplitude  y  is 
represented  in  terms  of  power  or  voltage. 

The  purpose  of  this  appendix  is  to  derive  an  expression  for  the 
expected  value  of  the  sample  variance  of  the  frequency-averages  of  the 
clutter  returns  from  the  N  resolution  cells  contained  in  the  range-azimuth 
stencil  illustrated  in  Fig.  A.l,  as  a  fui.tion  of  the  number  of  frequencies 
averaged,  stencil  size,  and  the  statistics  of  the  intrinsic  spatial 
probability  distribution  of  clutter  amplitude. 

Define  the  following  terms: 

The  sample  average,  xj_,  over  frequency,  of  the  returns  from 
cell  i: 
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yij 


(A. 3) 
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The  sample  variance,  V,  over  the  stencil  of  the  frequency-averaged 
returns  from  the  individual  cells: 
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(A- 5) 


Furthermore,  using  (A-2)  it  follows  that  for  a  given  cell  i: 

EZityij]  -  z±  (A. 6a) 

Therefore,  it  also  follows  from  (A. 3)  that: 

Ezitxil  *  zi  (A. 6b) 

where  the  notation  EZl  signifies  conditional  expectation 
given  Zj[. 
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Also  define: 


l(yirzi^  = 


V'ij 


°y/z1 


(A. 7)  represents  the  variance  of  the  frequency-diverse  returns  from  a  given 
resolution  cell,  the  returns  from  which  have  a  mean  value  z^.  It  should 

2 

be  noted  that,  in  general,  may  be  a  function  of  z^  and  is  therefore 

not  necessarily  constant. 

Finally  the  statistics  of  the  underlying  spatial  clutter  distribu¬ 
tion  are  defined  by: 

(A. 8a) 


„2  -  El*,) 


°z  ’  El(V“z)Z| 


(A. 8b) 


Given  the  preceding  definitions,  we  will  now  derive  an  expression 
for  the  expected  value  of  the  sample  variance  of  the  frequency  averages 
E [V] ,  equation  (A  ,5),  in  terms  of  the  number  of  cells,  N,  the  number  of 
frequencies,  J,  the  variance  of  the  frequency-diverse  returns, 


2 

°y/ zi * 


and  the  spatial  clutter  amplitude  statistics,  pz  and  az^» 


Taking  the  expected  value  of  equation  (A. 5),  one  can  write: 
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(A.  9) 
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Expand  Che  inner  term  of  the  above  summation  as  follows: 

\  [(>Ii '  *>2]  ■E2l[<xi '  zi> '  zi> 2] 

=  Ez  £  (x1  -  zi)2  -  2  (xA  -  z1)(X  -  z +  (X  -  z i)2J 


(A. 10) 


Now  expand  the  expectations  of  each  of  the  three  terms  in  (A.  10) 
individually  as  follows: 

For  the  first  term,  using  the  definition  of  x^  from  equation 


■ixi  ■  *i)2]- 3*  e  J  j,  (yii'2i)] 


Expanding  the  square  of  the  summation  and  talcing  the  expectation  inside  the 
resulting  double  summation  gives: 

\[<xi '  zi,2]'7  j, \  [(yu  ’  zi)<yit  ’  zi)] 

Since  the  frequency  samples  are  assumed  to  be  independent  it  follows  that: 

r  zi)(>,ik  ■  zi>]  ■  \(yij  -  zi>  -  zi>  -  °  f°r  j*k 

E2l[<yij  ■  zi)<yik  -  v]  ■  \[(yu  -  zi)Z]  for  j_k 
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Since  j“k  exactly  J  times  in  the  above  double  summation,  it  follows  that: 


E  (x  -  2  )2  =  I,  J  a.  2 
Zi  1  1 
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(A* 11) 


For  the  second  term  in  (A. 10) 
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Now  since  the  resolution  cells  are  independent,  it  follows  that: 

Ez  ^  ^Xi  ~  Zi^XA  ”  Zi  “  Ez  ^xi  “  zi^  Ez  for  1  = 


but  from  (A. 6b):  E  (x.  -  z .  )  **  0 
Zi  1  1 


Thus  Ez  (x±  -  z^x^  ~  z^)J  *  0  for  1  ”  i- 
a»d  [Ul  -  z1X«l  -  2t)]-  EzJ(«1  -  z4)2] 


for  i  *  £. 


Finally,  using  (A.  11) 


\  [<xi  ‘  zi,Cx*  '  zi>]  ■  T  °y2/zt 


for  i  ■  i. 


Substituting  this  back  into  the  expression  of  the  second  term: 
-2E 


^[<*1  -  z^a  -  v]-  -  IU-  »y/21 
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(A. 12) 


Finally,  for  the  third  term  in  (A. 10). 
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,  Ez  !*"  (x£  ‘  zi)2l 

'i  N  iLo,=l  1  1  J 


As  before,  expanding  the  square  of  the  summation,  and  taking  the 
expectation  inside  the  resulting  double-summation: 


Ez.  (X  -  V2  =  77  = 
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(A. 13) 


The  expectation  inside  the  double  summation  is  evaluated  individu¬ 
ally  for  each  of  the  following  four  cases: 


Case  A:  i  4  k  4  £  and 

Case  B:  i  *  k  4  £  and  i  =  £  f  k 

Case  C:  i  ■  k  *  £ 

Case  D:  i  4  k  *  £ 


For  Case  A:  i  4  k  4  £ 

EZi[(X£  •  Zl)(Xk  ~  *i>]“  (Z£  "  Zi)(Zk  '  *i> 


This  tase  occurs  N^-n-2(N-1)  times  in  the  summation  of  equation 


(A. 13) . 

For  Case  B:  i  *  k  4  £  and  i  4  £  =  k 

Ezt[(\  *  Zi)(xk  ‘  *!>]“  EzJ(xi  "  Zi)(xk  ‘  Zi>] 

*  Ezt(xi  -  zi)E(xk  ■  zi}- 


since  x^  and  x^  are  independent 


But  E(x£  -  z±)  ■  0 

Thus  EZi  (x^  -  z^Xx^  -  zj)  -  0 

this  case  occurs  2 ( N— 1 )  times  in  the  summation  of  equation  (A. 13). 

For  case  C:  i  *  k  *  i 

\  [  (x«  '  2l)(xk  -  21>]  '  \[  (xl  '  xl>2]-  7 

where  the  last  step  follows  from  (A*ll). 

This  case  occurs  only  once  in  the  summation  of  equation  (A-13). 

For  case  D:  14  (k=£) 

\  [<xt '  V  <xk  -  v] '  E2J<xi  -  xi>2]  w* 

Since  the  indicated  expectation  is  over  frequency,  expand  the  above 
as: 


Using  (A. 11)  and  (A. 6b),  gives: 

\K-  xi)2]-7  v*/  (xi’  v2 

This  case  occurs  (N-l)  times  in  the  summation  of  equation  (A. 13). 
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Now,  substituting  the  results  from  each  of  the  above  four  cases  into 


(A. 13)  gives: 
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Now  substituting  (A  .14),  (A  .12)  and  (A  .11)  into  (A  .10)  gives: 

vk  -  M>2]  - 1  -  -  °- 


v«t  W  °y/Zl  77  v.4  77  y/z! 
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k  £  4  N  £4i  * 


k*£*i 

Collecting  terms  and  substituting  into  (A. 9)  gives: 
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Expanding  the  expectations  In  the  last  two  summations  of  (A  .15)  gives: 

E  <v*i)(*r*i)  aEl[(2A)"(zA)ll(ziV‘(zA))l 
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Since  Is  uncorrelated  with  z^,  z ^  for  i*k,£  it  follows  that: 
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(A. 16) 


Now  recall  that  the  number  of  terms  in  the  double  sum  over  k,  £  for  k=£*i  in 
(A. 15)  is  n2-N-2(N-1)  •  (N-l)(N-2),  and  that  the  number  of  terms  in  the 
last  sum  for  £=*i  is  (N-l).  Therefore  when  the  result  of  (A  .16)  is  substi¬ 
tuted  into  the  summation  of  (A-15)  the  following  expression  is  obtained: 

E[V]  -  ^  L  E  1 0J  j+i.  I  [  (N-l)  (N-2)  a  2  +  2(N-l)o  *] 

N  J  i-1  y/Zi  NJ  i-1  Z  Z 

Finally,  performing  the  indicated  summation  gives: 
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Thus  the  expected  value  of  the  sample  variance  of  the  cell  frequency-aver¬ 
ages  over  the  stencil,  defined  by  equation  (A. 5)  consists  of  two  components 
as  indicated  in  equation  (A.17).  The  first  component  is  due  to  the  variance 
of  the  returns  from  any  given  cell  over  frequency,  and  this  component 
decreases  inversely  as  the  number,  J,  of  the  frequencies  averaged.  The 
second  component  is  the  underlying  variance,  over  space,  of  the  frequency- 
averaged  clutter  returns  from  the  different  cells  in  the  stencil.  This 
second  component  depends  upon  the  type  of  clutter  -  being  low  for  "smooth" 
clutter  such  as  flat,  mown  fields  and  high  for  "rough"  clutter  Buch  as  trees 
and  mixtures  of  two  or  more  different  clutter  types.  Note  that  this  second 
component  is  independent  of  the  number  of  frequencies  averaged,  and  effec¬ 
tively  places  a  lower  bound  on  the  expected  value  of  the  sample  variance,  V. 

Hence  the  effectiveness  of  frequency  averaging  depends  upon  the  rel¬ 
ative  magnitudes  of  these  two  components.  The  larger  the  variance  of  the 
intrinsic  spatial  amplitude  distribution  of  the  clutter  is  in  relation  to 
that  of  the  returns  over  frequency,  the  less  effective  frequency-averaging 
will  be  in  reducing  E[V], 

The  preceding  results  were  derived  without  making  any  assumptions 
about  the  specific  form  of  the  underlying  distributions  p(y/z)  and  p(z). 
However  in  order  to  apply  equation  (A.17)  to  real  data,  an  understanding  of 
the  actual  distributions  involved  is  helpful. 
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Specifically,  it  has  been  verified  by  Stovall^]  that  homogeneous 
clutter  (i.e.,  clutter  of  a  given  type  -  grass,  trees,  scrub,  etc.  -  in  con¬ 
trast  to  mixed  clutter)  can  be  very  well  represented  by  a  log-normal  distri¬ 
bution  whose  mean  and  variance  are  determined  by  the  specific  clutter  type 
in  question,  and  furthermore  that  the  distribution  of  the  returns  from  a 
given  resolution  cell  over  frequency  is  approximately  Rayleigh  or  exponen¬ 
tial,  depending  upon  whether  the  return  is  expressed  in  terms  of  voltage  or 
power  respectively.  Furthermore,  in  order  to  provide  the  maximum  possible 
dynamic  range,  logarithmic  detectors  are  often  used  for  amplitude  CFAR  sys¬ 
tems  of  the  type  being  analyzed.  As  a  consequence  of  the  logarithmic  detec¬ 
tion  process,  the  distribution  of  the  log-normally  distributed  clutter 
returns  will  appear  normal,  with  a  mean  and  variance  that  are  directly 
expressible  in  dB  (since  the  output  of  a  log  detector  is  linear  in  dB). 

It  is  shown  in  Appendix  C  that  the  variance  of  the  output  of  a  log 
detector  in  response  to  an  input  signal  that  is  Rayleigh  or  exponentially 
distributed  is  a  constant  regardless  of  the  mean  value  of  its  input.  As  a 
result,  when  the  frequency  diverse  returns  yjj  from  a  giver,  resolution 
cell  (cell  i)  are  obtained  as  outputs  of  a  logarithmic  detector,  the  condi- 

2 

tional  variance  Qy/Z  of  the  log  detector  outputs  will  be  constant , 
independent  of  their  mean  value  z^. 
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APPENDIX  B 


Normality  of  the  Probability  Distribution  of  the  Frequency  Averages  of 
Logarithmically  Detected  Clutter  Returns 


The  analysis  of  clutter  false  alarm  probability  in  section  3  was 
based  upon  the  assumption  that  when  the  returns  from  log-normally  distrib¬ 
uted  clutter  are  logarithmically  detected,  the  frequency  averages  of  the 
detector  outputs  will  be  normally  distributed.  The  purpose  of  this  appendix 
is  to  justify  this  assumption. 

It  has  been  shown  [Stovalll  that  the  amplitude  of  clutter  returns 
can  be  characterized  by  a  two-dimensional  probability  distribution  over  both 
radar  frequency  and  space  (resolution  cells)  in  the  following  way: 

py.«y(y-V  ‘  py/My(y/V  p«y(V 


where  y  is  the  amplitude  of  the  clutter  return  from  any  resolution  cell  at 
any  frequency  and  My  is  the  mean  value  of  the  clutter  return  from  a  given 
resolution  cell,  averaged  over  frequency. 

Thus  p(y/My)  is  the  conditional  distribution,  over  frequency,  of 
the  amplitude  of  clutter  returns  from  any  given  resolution  cell,  for  which 
the  mean  amplitude,  over  frequency,  is  My.  Hence,  My  represents  the 
intrinsic  amplitude  of  the  clutter  within  the  given  resolution  cell.  It  has 
been  postulated  and  experimentally  verified  that  p(y/My)  is  Rayleigh  dis¬ 
tributed,  given  the  specified  mean  value.  The  distribution  p(My),  on  the 
other  hand,  is  the  distribution  of  the  intrinsic  clutter  amplitude  over 
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space.  Hence  it  represents  the  distribution  over  space  of  the  frequency 
averages  of  clutter  returns  from  any  given  resolution  cell.  It  has  also 
been  verified  that  p^CMy)  is  log-normally  distributed  for  a  wide 
variety  of  clutter  types. 

This  two-dimensional  separation  of  the  distribution  of  the  clutter 
returns  over  both  frequency  and  space  is  valid  regardless  of  any  intermedi¬ 
ate  functional  transformation  of  the  returns.  Hence  if  the  y's  represent 
the  output  of  a  logarithmic  detector,  the  previous  argued  separability  of 
the  two-dimensional  distribution  still  holds,  except  that  now  the  form  of 
Py/My(y/My)  will  be  logged-Rayleigh  while  that  of  PMy(My)»  will 
be  normal.  However  since  the  x^'s  are  suns  of  independent  r.v.  's  (the 
log-detector  outputs  yjj)  each  of  which  are  distributed  in  accordance  with 
Py/My(y/My),  it  follows  from  the  central  limit  theorem  that  for  a 
sufficiently  large  number  of  frequencies  averaged,  the  xj/s  will  all  be 
approximately  normally  distributed.  This  is  particularly  true  if  the 
yjj's  are  logged-Rayleigh  distributed  as  postulated,  since  that  distribu¬ 
tion  has  a  shape  similar  to  that  of  a  normal  distribution  of  like  mean  and 
variance.  Because  of  the  separability  argument  the  distribution  of  the 
Xj_'s  may  again  be  written: 

p(x,Mx)  -  Px/i^)  PMjjCMx) 

Furthermore,  since  the  mean  of  the  sum  is  equal  to  the  sum  of  the  means,  it 
also  follows  the  M*  ■  My.  Thus  the  distribution  of  p^CM*)  will 
be  the  same  as  PMy(My)  -  i.e.,  normal  for  the  case  of  logarithmic 
detection  of  log-normally  distributed  clutter  returns,  while  that  of 
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Px/b^Cx/Mx)  will  be  approximately  normal  by  virtue  of  the  averaging 
process  used  to  form  the  x^’s  and  the  central-limit  theorem.  Hence  the 
density  p(x,Mx)  is  the  product  of  two  normal  distributions  and  can  be 
shown  to  be  jointly  normal.  The  r.v.  ’s  x  and  M*  are  of  course  not 
independent.  The  distribution  of  the  x^  for  all  M*  is  then  given  by  the 
marginal  density  of  pCx.Mjj)  as  follows: 

P*(x)  =  /  P(x’Mx>  ^X  “  i  Px/M  (x/V  Pm  (Mx>  dMx 

X  X 

—oo  —oo 

Since  p(x,Mx)  is  jointly  normal,  the  density  px(x)  will  also  be  normal. 
This  completes  the  verification  of  the  normality  of  the  distribution  of  the 
xj  when  log-normal  clutter  is  logarithmically  detected. 
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APPENDIX  C 


Variance  of  the  Output  of  a  Logarithmic  Detector  in  Response  to  a 
Rayleign  Distributed  Noise  Input 


Let  v  be  the  magnitude  of  the  complex  signal  voltage  at  the  input  to 
a  logarithmic  detector.  If  v  is  a  Rayleigh  distributed  random  variable, 
then  its  probability  density  function  is  given  by: 


fv(v) 


2 o, 


(C.l) 


where  ax^  is  the  variance  of  the  real  and  imaginary  parts  of  the  underlying 
complex  Gaussian  random  process.  For  the  present  purposes  it  is  more  con¬ 
venient  to  work  with  the  instantaneous  power  corresponding  to  this  input 
voltage,  rather  than  with  the  voltage  itself.  Hence,  define: 

q  -  V2  (C.2) 

The  distribution  of  q  is  easily  found  to  be  [C-l] : 

e"q/uq 


fq(q) 


Vq 


(C.3) 


where  Pq  *  2a the  mean  value  of  v2  -  the  mean  input  power.  The  mean 

value  of  the  input  voltage  magnitude,  v,  is  related  to  pq  by:  v  “^"^Pq 

The  output  of  a  logarithmic  detector  can  be  calibrated  to  be  linear  in  dB 

as  follows:  let  x  represent  the  voltage  at  the  output  of  the  logarithmic 

detector,  then: 


x  ■  K  log.nv,  volts 


(C.4) 


where  K  is  the  log  detector  proportionality  constant.  However,  the  decibel 
level  of  the  input  signal,  relative  to  an  arbitrary  1  volt  level,  is 
defined  by: 

y  -  20  log1Qv,  dB  (C.5) 

_  20 
Thus  y  “  x 

In  other  words,  multiplying  the  log-detector  output  voltage,  x,  by  the 
calibration  constant  (20/K)  yields  the  signal  level  in  dB. 

Substituting  equation  C.2  into  C.5  gives: 

y  ■  io  log 10  q 

Furthermore,  since  log^g  q  =»  0.4343  Ln  q  it  follows  that: 

y  ■  4.343  Ln  q  dB  (C.6) 

We  wish  to  determine  the  variance  of  y  for  any  Rayleigh  distributed  input 
signal  of  arbitrary  magnitude  (mean  power  uq,  or  mean  voltage  /pq  ) 

By  definition: 

ay  -  E[y2]  -  e2[7]  (C.7) 

Using  the  density  function  of  q,  given  by  equation  C.3,  and  the  definition 
of  y  given  by  C.6  it  follows  that: 


®  e-cl/vq 

E[y]  =  4.343  f  (Ln  q)  -  dq 

o  ^  q 

®  o  e-q/^9 

E[y2]  “  (4.343)2  J  (Ln  q;  - dq 

o  ^q 

The  integrals  in  C.8a  and  C.8b  are  tabulated  in  Ref.  C.2: 

“  e"q/wq 

/  Ln  q - dq  =  -0.57721  +  Ln  yq 

o  wq 

00  9  Q-q/Pq  o  _2 

f  (Ln  q) £  - - dq  =  (-0.57721  +  Ln  pq)  + 


Thus  E[y]  -  4.343  (-0.57721  +  Ln  pq)  (C.9a) 

2 

E[y2]  =  (4.343) ^  [(-0.57721  +  Ln  uq)2  +  \  ]  (C.9b) 

2 

Substituting  these  results  into  the  definition  of  a y  equation  C.7,  gives: 

Oy2  =  (4.343)Z  =  31.025  dB* 


or  Oy  =*  5.57  dB 

Note  that  this  result  is  independent  of  the  mean  value  of  the  Rayleigh 
distributed  input  signal. 
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